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Abstract 



An anomaly-free operator corresponding to the Wheeler-DeWitt constraint 
of Lorentzian, four-dimensional, canonical, non-perturbative vacuum gravity 

■ is constructed in the continuum. This operator is entirely free of factor order- 
ing singularities and can be defined in symmetric and non-symmetric form. 

We work in the real connection representation and obtain a well-defined 

■ quantum theory. We compute the complete solution to the Quantum Einstein 
Equations for the non-symmetric version of the operator and a physical inner 

q \ product thereon. 

\Q • The action of the Wheeler-DeWitt constraint on spin-network states is 

by annihilating, creating and rerouting the quanta of angular momentum 
\ associated with the edges of the underlying graph while the ADM-energy is 

essentially diagonalized by the spin- network states. We argue that the spin- 
Q" 1 network representation is the "non-linear Fock representation" of quantum 

', gravity, thus justifying the term "Quantum Spin Dynamics (QSD)". 



1 Introduction 

Attempts at defining an operator corresponding to the Wheeler-DeWitt constraint 
[0] of Lorentzian, four-dimensional canonical vacuum gravity have been made first 
in the framework of the metric (or ADM) variables (see, for instance, ||). The 
formulation of the theory seemed hopelessly difficult because of the complicated, 
non-polynomial algebraic form of the Wheeler-DeWitt constraint and it was there- 
fore thought to be mandatory to first obtain a polynomial formulation of the theory 
to complete the programme. 

The celebrated observation due to Ashtekar [|[] is that the canonical constraints 
of general relativity can indeed be cast into polynomial form if one performs a certain 
complex canonical transformation on the gravitational phase space. This opened, 
for the first time, the hope that one can actually rigorously define the quantum 
Hamiltonian constraint (or Wheeler-DeWitt equation). 
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The major roadblock after this was the difficult reality structure of the theory so ob- 
tained : general relativity, when written in Ashtekar's variables, is a dynamical the- 
ory of complex- valued connections for the non-compact gauge group SL(2, (D), how- 
ever, standard mathematical constructions and techniques usually used for Yang- 
Mills theory apply only if the gauge group is compact. 

A solution to this problem has been recently proposed in form of a phase space 
Wick rotation transform [El] (see also ||) which enables one to formulate general 
relativity as a SU(2) gauge theory while keeping the polynomial algebraic form 
of Wheeler-DeWitt constraint and incorporating the correct reality conditions into 
the quantum theory. This theory can be, not surprisingly, recognized as Euclidean 
four-dimensional gravity as written in terms of real valued Ashtekar variables @. 
This then opens access to the recently developed, powerful calculus on the space of 
(generalized) connections modulo gauge transformations for compact gauge groups 
[§, |9|. [U], [11], This calculus provides a rigorous kinematical framework by means 
of which constraint operators can be regularized in a well-defined and unambiguous 
manner. In particular, this framework has already been successfully employed to 
arrive at the general solution of the Gauss and Diffeomorphism constraints fl3| . 

In order to complete the programme one still needed to construct a rigorously 
defined operator corresponding to the classical generator of the transform. This 
seems a hopeless problem to solve as the classical generator is a non-polynomial, 
not even analytic function of the phase space variables. 

There is, however, an even more severe problem : the Hamiltonian constraint 
of Euclidean or Lorentzian gravity as defined by Ashtekar is a density of weight 
two. Namely, in order to obtain a polynomial form of those constraints one needs 



to rescale the constraint functional and to absorb a factor of l/ydet(g) into the 
Lagrange multiplier (the lapse function; here q = (q a b) is the intrinsic metric of 
an initial data hyper surf ace). On general grounds, an operator corresponding to a 
classical function on the phase space with density weight different from one needs 
to be renormalized. While this poses no, a priori, problems for, say, Yang-Mills 
theories in a classical background geometry, such a procedure is unacceptable for 
quantum gravity since a renormalization unavoidably introduces a length scale and 
thus breaks diffeomorphism invariance. 

A solution to this problem was first suggested in jL4] and is currently reconsidered 
in : to produce an operator of density weight one one takes the square root of 
the rescaled Euclidean Hamiltonian constraint. While classically this does not alter 
the theory and while this seems to be a necessary step to do in quantum theory 
in order to keep diffeomorphism invariance, there remain problems that have to do 
with taking the square root of an infinite number of non-self adjoint, non-positive, 
non-commuting operators. 

In the present article we suggest a technique to solve both problems, namely the 
complicated reality structure of the Lorentzian theory as well as the problem con- 
nected with the density weight, in one stroke : what was thought to create problems 
in the quantization process is precisely the reason for why we are able to find a finite 



operator : the factor l/\/det(g) is needed. 



• We show that it is indeed possible to construct a finite, (symmetric) operator 
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corresponding to the original, non-rescaled, Lorentzian Wheeler-DeWitt con- 
straint whose quantum constraint algebra is non-anomalous. Since the original 
Wheeler-DeWitt constraint has density weight one, no renormalization is nec- 
essary. 

Amazingly, the resulting operator is not messy at all and the problem of finding 
exact solutions to the quantum constraint is conceivable. 

We always work with real-valued Ashtekar variables, the reality structure of 
the theory is very simple, the complex- valued Ashtekar variables are never 
introduced. 

In the construction of the Wheeler-DeWitt operator the unrescaled Euclidean 
Hamiltonian operator as well as the generator of the Wick transform arise 
in a natural way so that both operators are also constructed rigorously as a 
side result. This is important since the Wick rotation transform simplifies the 
problem of finding solutions to the quantum constraint. 

Using the same technique one can give rigorous meaning to a whole bunch of 
other operators in a representation where the metric is not diagonal, including 



but not exhausting a) the operator corresponding to the length of a curve [16 
b) the generators of the Poincare group for asymptotically flat topologies |f[7 
and c) matter contributions to the Hamiltonian constraint E8|. It might be 
that it is in this sense that quantum gravity arises as the "natural regulator 
of the matter field theories" . By this we mean the following : in the canonical 
quantization programme a natural regularization procedure is point splitting. 
With the exception of spinorial matter, all matter Hamiltonians of, say, the 
standard model, are quadratic in the momenta which means that they display 
a density weight of two when neglecting the gravitational interaction while they 
have density weight one when taking gravity into account. When removing 
the regulator the density weight shows up in the form of a product of delta 
distributions evaluated at the same point which is singular. On the other 
hand, with gravity the singularity is removed, in the limit one arrives at a 
well-defined operator-valued distribution. 

The article is organized as follows : 

After fixing the notation and explaining the main idea we construct first the 
unrescaled Euclidean Hamiltonian constraint operator. We motivate our choices 
involved in the regularization step. The freedom in our choices is severely restricted 
by the requirement that the resulting operator be diffeomorphism-covariantly defined 
and non-anomalous. We provide a solution to both requirements. If one adds the 
requirement that the operator be at least symmetric (preferrably, it should possess a 
self-adjoint extension) then the regularization involves an additional structure. We 
will stick with a non-symmetric operator in the main text and provide a symmetric 
operator in jl| which is non-anomalous and diffeomorphism-covariantly defined as 
well. 

Next we construct operators corresponding to the generator of the Wick rotation 
transform and finally the Lorentzian Wheeler-DeWitt operator. 
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For the non-symmetric operator we are able to find the complete kernel of the 
Wheeler-DeWitt constraint operator in Q. The physical Hilbert space turns out to 
be the one already given in [ p~3|j . 

For the symmetric operator on the other hand we do not have the complete so- 
lution yet, although solutions can be easily computed by a case by case analysis. In 
principle, since, expectedly, on diffeomorphism invariant states the constraint alge- 
bra is Abelian we are able to find solutions to the constraint as well as constructing 
a physical inner product by the group averaging method [19, ^(J. We also comment 
on finding observables along the line of argument in |0| . A lot of the mathematical 
problems of quantum gravity can be solved, on so-called cylindrical subspaces, with 
well-known Hilbert space techniques familiar from quantum mechanics. 

We have complete control over the space of solutions to both versions of the 
constraint and the intuitive picture that arises is the following : The Hamiltonian 
constraint acts by annihilating, creating and re-routing the quanta of angular mo- 
mentum (with which the graphs of so-called spin- network [21, 22, 23 states are 
"coloured") in units of ±h, ih/2, 0. On the other hand, linear combinations of such 
states diagonalize the ADM energy operator very much in the same way as Fock 
states diagonalize the Maxwell Hamiltonian, the role of the occupation number being 
played by the spins of the spin-network state. Thus, the spin-network representation 
is the "non-linear Fock representation for quantum gravity" and this motivates to 
call the quantum theory we obtain "Quantum Spin Dynamics (QSD)" in analogy 
with QED or QCD. 



2 Notation and the main idea 

Let the triad on the spacelike, smooth, hypersurface S be denoted by e l a , where 
a, b, c, ... are tensorial and k, ... are 577(2) indices. The relation with the intrinsic 
metric is given by q ab = e, l a e b^ij- ^ follows that det(g) := det((g a &)) = [det((ejj)] 2 > 
0. The densitized triad is then defined by Ef := det((e£))ef where e" is the inverse 
of e\. We also need the field K l a = eK ab e b , e = sgn(det((e£))), where K ab is the 
extrinsic curvature of E. It turns out that the pair (K l a , Ef) is a canonical one, that 
is, these variables obey canonical brackets {K^(x), E h Ay)} = k8^ (x,y)S b Sj where k 
is Newton's constant. 

Let the spin-connection (which annihilates the triad) be denoted by r*. Then one 
can show that (A l a := + K l a , Ef) is a canonical paiif] on the phase space of 
Lorentzian gravity subject to the SU(2) Gauss constraint, the diffeomorphism con- 
straint and the Wheeler-DeWitt constraint (neglecting a term proportional to the 
Gauss constraint) 

H := ^/det(q)[K ab K ab - {K a a f - R] = 1 ti{{F ab - 2R ab )[E a , E b ]) (2.1) 

v /det(g) 

where F ab and R ab respectively are the curvatures of the SU (2) connection A l a and 
the triad e\ respectively. 

1 If we had chosen Ef = y/dct(q)ef instead of Ef = det(e)e" then this pair is not canonical, a 
fact often overlooked in the literature 0. Roughly speaking, we would spoil the integrability of 
the spin connection. 
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What has been gained by reformulating canonical gravity as a dynamical theory 
of 577(2) connections is the following : if, as we do in the sequel, one makes the 
assumption that there exists a phase for quantum gravity in which the excitations 
of the gravitational field can be probed by loops rather than, say, test functions of 
rapid decrease, then one has access to a powerful calculus on the space of (gener- 
alized) connections modulo gauge transformations A/Q and, in particular, there is 
a natural choice of a diffeomorphism invariant, faithful measure //o thereon which 
equips us with a Hilbert space TC := L 2 (A/Q, dfi ), appropriate for a representation 
in which A is diagonal. Moreover, Gauss and diffeomorphism constraints can be 



solved (see [10] and references therein for an introduction to these concepts). 
The remaining step then is to give a rigorously defined quantum operator corre- 
sponding to the Wheeler-DeWitt constraint and to project the scalar product on its 
kernel. 

Let us explain elements of the underlying kinematical framework [Of : 



We begin by explaining the notion of a "cylindrical function" on A/Q. In brief 
terms, gauge invariant cylindrical functions on the space of (generalized) 577(2) con- 
nections are just finite linear combinations of traces of the holonomy around closed 
analytic^ loops in E. Each such function therefore may equally well be labelled by 
the closed, piecewise analytic graph 7 consisting of the union of all loops involved in 
that linear combination. Such a graph consists of a finite number of edges ei, .., e n 
and vertices V\, ..,v m . So, a function cylindrical with respect to a graph 7 typically 
looks like f(A) = / 7 (/i ei (A), ..,h en (A)), where h e (A) is the holonomy along e for 
the connection A and f 1 is a complex- valued function on SU(2) n such that f(A) 
is gauge- invariant. The functions cylindrical with respect to a graph that are n- 
times differentiable with respect to the standard different iable structure on SU(2) m 
for some m will be denoted by CjV^{A/Q) and Cy\ n (A/Q) := U^Cyl^A/Q). One 
and the same cylindrical function can be represented on different graphs leading to 
cylindrically equivalent representants of that function. It is understood in the above 
union that such kind of functions are identified. 

The measure /io referred to above is entirely characterized by its cylindrical projec- 
tions defined by 

fdno(A)f(A)= fdn „(A)f^({h ei (A)}) = [ d n {i H (g u .., g n )f y (g u ••, g n ) ■ 
JA/g JA/g Jsu{2)" 



An orthonormal basis on 7i is given by the so-called spin-network states [El], |22], ^] 
: given a graph 7, "colour" each of its edges e with a non-trivial irreducible repre- 
sentation 7Tj e of 577(2), that is, j e is the spin associated with e. With each vertex 
we associate a contraction matrix c v which contracts all the matrices iTj e (h e ) for e 
incident at v in a gauge-invariant way. In this paper we will denote a spin-network 
state by j s if we wish to stress the dependence on 7, j = (j e ), c— (c v ) where the 

vectors j, c have indices corresponding to the edges and vertices of 7 respectively. 
Consider the set of smooth cylindrical functions $ := Cyl°° (A/Q) which can be 
shown to be dense in 7i. By a distribution tp 6 $' on $ we mean a generalized 
function on A/Q such that for any (f> G $ the number i/)[<f>] := j-jjg dn ip(f) < 00 is 

2 The extension of the framework to the smooth category is possible but not entirely straight- 
forward |E4) . For simplicity we restrict ourselves to the analytic category. 
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finite. It turns out that the solutions of the diffeomorphism constraint are elements 
of $' @. 

Suppose we have a quantization H(N) of the Hamiltonian constraint operator, that 
is, it is densely defined on 7i and its classical limit reduces to H(N), where N is 
the lapse function. Then its adjoint H(Ny also has the same classical limit H(N) 
(a reordering of terms gives only higher orders in h) because H(N) is real- valued. 
Therefore, if H(N) is not self-adjoint, then it is an option of whether we impose 
H(N)ip = or H^Nyip = 0. In both cases, the solution ip is typically not an L 2 func- 
tion any longer but a distribution, that is, in our case an element of Thus, given 
the framework of generalized eigenvectors |13| and the associated triple $ C Ti C $' 
we choose to find the generalized eigenvectors with eigenvalue corresponding to 
the kernel of the (self-adjoint) Hamiltonian constraint operator H(N)^ as follows : 
Let tp G be a distribution. We say that ip is in the kernel of H(Ny whenever 
(H(Nyijj)(f) := ip(H(N)f) = for each lapse and each cylindrical function / in 
the (dense) domain of H(N). Note that we cannot require that ip is diffeomorphism 
invariant if we impose the condition H(Nyip = 0. This is because the Hamilto- 
nian constraint does not leave the subspace of corresponding to diffeomorphism 
invariant elements, invariant so that one would be forced to solve the Hamiltonian 
constraint before the diffeomorphism constraint. On the other hand, we will see that 
if ip is diffeomorphism invariant and H{N) is diffeomorphism covariantly defined, 
then solving ip(H(N)f) = is meaningful. 

After this preparation we are now ready to explain the main idea of our approach. 
Suppose that we can give meaning, in a representation in which A is diagonal, to 
two operators corresponding to 

1) The total volume of £ given by[] 

V := V(E) := J d 3 xyj\det{q)\ and (2.2) 

2) the integrated trace of the (densitized) extrinsic curvature of £ 

K := J d 3 x^Jdet(q)K ab q ab = J d 3 xK l a E? . (2.3) 

This means that their quantizations V, K are densely defined on suitable subspaces 
of cylindrical functions (in case that we wish to obtain a symmetric operator we 
will also require that they are self-adjoint on 7i). The motivation for these two 
assumptions comes from considering the following two key identities 

1 . ' J = e abc eei(x) = 2e abc -—-- = 2e abc {^, V} and (2.4) 
yjdet(q) 6E i( x ) K 

Ki = -r— = { — ,K\. (2.5) 



3 In case that E is not compact, as we will need only the variation of V in the sequel, (2.2) 
is to be understood in the following way : consider any nested one-parameter family of compact 
manifolds T, r C E r <, V0<r<r'<oo such that lim^oo E r = E. Then 5V(Y<) := linv^oo 5V(Y, r ) 
and this is well-defined. 
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The last equality relies on the observation that {T^, K} = and it is this identity 
underlying the ideas developed in ||. The importance of these identities becomes 
clear when we get rid of the complicated curvature term R ab involved in ( |2.1| ) in 
favour of K. We have 

H + H E = - 7 =£=tr([K a , K b \ [E\ E b }) = -— L=tr([ {^, K}, A K}} [E\ E b }) 
v /det(g) yjdet(q) « « 

= ±e ab %r([{A a , K}, {A b , K}}{A C , V}) = ^e abc tr({A a , K}{A b , K}{A C , V}) (2.6) 
where we have introduced the (unrescaled) Euclidean Hamiltonian constraint 

H E := , 1 tr(F ab [E a , E b ]) = -e abc tr(F ab {A c , V}) . (2.7) 
/det(g) k 

So what we have achieved is to hide the non-polynomiality of the theory as de- 



termined by l/ydet(g) in a Poisson bracket. Classically this is not helpful at all 
(except, possibly, for Hamilton- Jacob i methods or semi-classical approximations), 
however, we will show that in the quantum theory it is of advantage. Namely, it 
is now clear where we are driving at : in any regularization of the Wheeler-DeWitt 
constraint operator we will have to approximate the connection A a and the curva- 
ture F ab respectively by cylindrical functions given by the holonomies h along some 
edges and closed loops respectively. Now, one obvious quantization of ( |2.6j ) would 
be to replace V, K by V,K and the Poisson brackets {.,.} by [.,.}/ (ih). It then 
follows, given our assumption, that this quantization of ( |2.6|) has a chance to result 
in a finite operator on cylindrical functions since no singular terms appear when 
computing the bracket and we would have managed to produce a densely defined 
operator. 

Is it then true that there exist quantizations of V, K meeting our assumption ? The 
answer is, surprisingly, in the affirmative : 

First, it is a fact that there has already been constructed a well-defined, self-adjoint 



operator V on Ti, corresponding to V p5| . p6| whose action on cylindrical functions 
is perfectly finite : (we follow |p6| ) 



Vf = 4 E V v f= Ul £ \— £ e{e u ej,e K )e ijk X\X^X\\ f 7 ( gi ,..,g n ) 

»eV(7) \ veV(-y) \ eiC\ejr\e K =v J 

(2.8) 

where e(ej, ej, ex) = sgn(det(e/(0), ej(0), eR-(O))). We have abbreviated gi = h ei (A) 
and Xi = X(gj) is the right invariant vector field on SU(2) (we have chosen orien- 
tations such that all edges are outgoing at v). V(j) is the set of vertices of 7 and 
l p := \fhn is the Planck length. The symbol V v is defined as follows : 



V f _ i ^\^Y,e I nejne K =v^I^J^K)e ij kX i I X 3 jX k K \ f y if V G V(j) 

vh ~ \ ' if«*V(7) ' 

that this, it is an operator which acts on the edges of any graph meeting at the 
point v in the way displayed in Q2.8| ). This demonstrates that V is a densely defined 
operator on thrice differentiable cylindrical functions. From this it follows already 
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that we have a chance of giving meaning to an operator corresponding to H E . 
Secondly, it is a well-known fact that K is, up to a multiplicative constant, just 
the time derivative of the total volume with respect to the integrated Hamiltonian 
constraint (which is a signature invariant statement) 

K = -{-, / d 3 xH E (x)} (2.9) 

which of course can also be verified immediately. So, if we (again) replace V, H E 
by their quantizations and Poisson brackets by (l/(ih) times) commutators we also 
find that we have a chance of giving meaning to an operator corresponding to K. 
This completes our explanation of the main idea. The rest of this paper is devoted 
to a precise construction of the operators sketched in this section. We do this in a 
series of three steps : 

Step A) We begin by giving meaning to an operator corresponding to the Euclidean 
Hamiltonian constraint (|2.7|) . The result will be a (self-adjoint) operator on 7i whose 
constraint algebra is anomaly-free. 

Step B) We quantize K along the lines sketched above using the known quantizations 
V ,H E . 

Step C) We quantize the Wheeler-DeWitt constraint using the known quantizations 
H E ,K and exploiting ( |2.6| ). We show that its constraint algebra is anomaly free 
and that the result is a (symmetric) densely defined operator on 7i as well. 
Recalling that C = (tt/2)K is the classical generator of the Wick rotation transform 
|| we naturally obtain its quantization, as well as that of the Euclidean Hamiltonian 
constraint, for free in our procedure. 



Quantization of the Euclidean Hamiltonian con- 
straint 



The method applied in Hl4| , jig] is to absorb the prefactor l/ydet(g) in (|2.7| ) into 
the lapse function and to give meaning to the operator corresponding to the square 
root of the trace. We will not do this. By employing the method described below we 
can avoid the complications that arise in connection with these two steps. Also, it 
would be less straightforward to define K with this rescaled form of the constraint 
since then K is not the time derivative of the volume any longer. 



3.1 Regular izat ion 

The regularization and all computations will be performed in an arbitrary but fixed 
standard frame for S as usual. The end result will be independent of that choice of 
frame. 

We start from the classical expression for the Euclidean constraint functional 

H E [N] = - [ d 3 xN(x)e abc tr(F ab {A c ,V}) (3.1) 

K JT, 

where N is the lapse function divided by k. 

We now triangulate E into elementary tetrahedra A where each of its edges are 
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analytic. For each tetrahedron we single out one of its vertices and call it v(A). 
Let Si (A), i = 1,2,3 be the three edges of A meeting at v(A). Let ctij(A) := 
Si (A) o ajj(A) o Sj(A) -1 , cKji = a~- , be the loop based at u(A) where is the 
obvious other edge of A connecting those endpoints of Sj, Sj which are distinct from 
v(A). Then it is easy to see that 



Hi[N] := -\N v e^\x{h^h sm {h- s l {A) , V}) (3.2) 



tends to the correct value 2 j A [Ntr(F A {A, V})] as we shrink A to the point 
v(A), N v := N(v(A)). Moreover, H A [N] is manifestly gauge-invariant since «y o 
Sk o SjT 1 is a "loop with a nose" . 
Let the triangulation be denoted by T. Then 

H§[N] = £ F|[iV] (3.3) 
Aer 



is an expression which has the correct limit (^.1[) as all A shrink to their basepoints 



(of course the number of tetrahedra filling any bounded subset of E grows to infinity 
under this process). 

As we have said before, we now simply replace V by V and the Poisson bracket by 
\j%h times the commutator and 

H*[N] := £ H*[N}, Hl[N] : = -2^^e^ fc tr(/ iQij(A) / isfc(A) [/ i ; fc 1 (A) , V}) =: N V H* 
Agt p 

(3.4) 

is an operator with the correct classical limit. 

It is obvious that the properties of the operator ( p.4|) are largely determined by the 
choice of triangulation T and therefore we devote the subsequent paragraph to a 
preliminary investigation of those properties of ( |3.4|) which hold for any choice of 
triangulation. These considerations will then motivate our choices. 



3.1.1 Motivation 

The first property of (|3.4j ) that we wish to prove is that its action on cylindrical 
functions is indeed finite no matter how fine the triangulation T is, provided that a 
certain criterion is satisfied which we derive now. 

In order to see this it is sufficient to consider the operator [/i Sfe (A), V]f where / is a 
cylindrical function with respect to some graph 7. The first case is that Sfc(A)fl7 = 0. 
Then Vfr U s fc (A)) = ^(7) U V(s k (A)) and it follows from O that [h Sk(A) , V}f = 

T, v ev(-r) h Sh (A)V v f — J2veV(yUs k (A)) Vvh Sh (A)f = — fJ2vev(s k (A)) V v h Sk (A) = since V 
annihilates any cylindrical function, whether gauge- invariant or not, whose under- 
lying graph is not at least three-valentQ. 

The next case is that Sfc(A)fl7 7^ but does not contain a vertex of 7. That is, the set 
V(7 Usfc(A)) — V(7) consists of vertices which are at most four-valent, however, the 
tangents of the edges incident at each of those vertices lie in a two dimensional vector 
space. It then follows that again [h Sk (A),V)f = - E ve y( 7 u Sfc (A))-v( 7 ) KA fc (A)/ = 

4 We say that a graph is n-valent if there are no more than n edges ingoing or outgoing at each 
vertex 
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because of the signature factor e(sj, Sj, Sk) involved in (|2.8|). Notice that this prop 
erty would no longer be true had we used the volume operator as defined in [p5 | 



according to the second reference in |215 |, that operator does not vanish for vertices 
which involve only edges with co-planar tangents. Therefore, had we used this op- 
erator, as we make the triangulation finer and finer we would get more and more 
contributions and thus the resulting continuum operator could not even be densely 
defined. 

It follows that (13.41) reduces to 



H§[N]f= E Hi[N}f= £ N v ^Hif=: £ N V H? f (3.5) 

AnV(7)^0 veV(y) veA veV{j) 

and we see that this expression is finite no matter how "fine" the triangulation is, 
provided the number of tetrahedra intersecting the vertices of 7 stays bounded as 
we go to finer and finer triangulations. These considerations motivate to construct 
a triangulation T{^f) assigned to a graph 7 which meets this criterion and we get an 
operator H E ,y This furnishes our preliminary analysis. 

3.1.2 Requirements for a triangulation adapted to a graph 

Certainly, there are an infinite number of possible assignments. We choose a partic- 
ular assignment guided by the following principles : 

a) Non-triviality : 

We could choose T(j) in such a way that there are no intersections with 7 at 
all, giving automatically a trivial result. This is inappropriate as the number of 
degrees of freedom should be genuinely reduced by the Hamiltonian constraint. 

b) Diffeomorphism-Covariance : 

Remember that we want to impose the constraint on a distribution if) as out- 
lined in section 2 and that the measure /xq is diffeomorphism invariant. This 
fact enables us to get rid of a huge amount of ambiguity arising in the assign- 
ment of a triangulation as follows : the classical Hamiltonian constraint is not 
diffeomorphism invariant but it is diffeomorphism covariant. If we could carry 
over this classical property to the quantum theory then diffeomorphic vectors 
would be mapped by H E into diffeomorphic vectors. Therefore, provided the 
state if) is diffeomorphism invariant, the number ip[H E f] would only depend on 
the diffeomorphism invariant properties of the triangulation assignment (this 



was first observed in ||14J| ). The assignment should therefore move with the 
graph 7 under diffeomorphisms of E. 

More precisely, we have the following^] : Let 7 be a graph and 0(7) its diffeo- 
morphic image for some smooth diffeomorphism G Diff(E) (the diffeomor- 
phism does not need to be analytic but it must keep the graph analytic). Then 
Hrpr^fi will depend on a graph 7 and likewise H^um) f<j>(y) w ih depend on a 

graph 0(7). Then the requirement is that 7,0(7) are diffeomorphic. 

That is, we want that there exists 0' G Diff(E) such that U{<f/)[H^/ sf^\ = 



5 This more precise formulation of the intuitive idea of diffeomorphism covariance was commu- 
nicated to the author by Jurek Lewandowski 



10 



H-TUM)f<t>('y) wnere U((fi')fry = f<p'(j) is a unitary representation of the diffeo- 
morphism group Diff(S) on 7i. Notice that we have been dealing with the 
Hamiltonian constraint at a point H E rather than with its smeared version 
H E [N] which is appropriate because we need to impose the constraint at every 
point of X. We will see that in our case the notion of a constraint operator at 
a point makes perfect sense. 

c) Cylindrical consistency : 

If 7' is bigger than 7 then T{^) and T(j') will in general differ from each 
other. However, if / is cylindrical with respect to 7 then the vectors H E ^,^f 

and H E ^f should be diffeomorphic to each other. That is, we have cylindrical 
consistency up to a diffeomorphism. The reason why we do not need to require 
exact cylindrical consistency is because the assignment of the triangulation is 
only defined up to a diffeomorphism if we care only about the evaluation of a 
diffeomorphism invariant state ip on the states H E ff\. 

d) Symmetry and Self-Adjointness : 

The classical Hamiltonian constraint is a real-valued function on the phase 
space of general relativity. It is therefore compatible with the principles of 
quantum theory to construct an operator corresponding to it which is self- 
adjoint or at least symmetric. While one is not forced to do so as the symmetric 
and non-symmetric operators have the same classical limit and as we are only 
interested in the point of the spectrum, rather than the full spectrum, there 
are practical reasons, among others the applicability of the group averaging 
method and the possibility of being able to get rid of a quantization ambiguity, 
which motivate to have the constraint in symmetric form. We will therefore 
propose two quantizations of the Hamiltonian constraint : in the main text we 
will stick with a non-symmetric operator as it turns out that it is technically 
much easier to handle and in |]J we will treat a symmetric version of the 
operator. 

e) Efficiency : 

The result of applying H E [N] to cylindrical functions will be a cylindrical 
function that depends on additional edges. We want to choose an assignment 
which introduces as less additional structure (edges) as possible. 

f) Naturality : 

The assignment should be uniform, that is, it treats all the edges of the graph 
7 incident at a vertex on equal footing. 

g) Anomaly- freeness : 

The assignment should be free of anomalies, that is, the constraint algebra 
should close, otherwise we are reducing the number of degrees of freedom too 
much and we do not obtain a quantization of the classical theory. 

h) Non-Emptyness : 

The assignment should not be such that the kernel of the constraint operator 

6 This observation again was first communicated to the author by Jurek Lewandowski 
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is empty. The space of solutions to the classical Einstein equations has a rich 
structure and so an empty kernel is not appropriate as it would not correspond 
to the non-empty classical reduced phase space. 

3.1.3 Choice of a triangulation adapted to a graph 

The number of choices meeting these requirements is certainly still infinite. The 
assignment T(j) we choose is as follows : we give in this subsection an assignment 
which is appropriate only for the non-symmetric regularization of the operator and 
will modify it at a later stage in order to make it appropriate for the symmetric 
version. 

0) Two-valent vertices : 

If a vertex v is two-valent, adjoin one more edge to 7 incident at v and not 
intersecting 7 in any other point such that its tangent at v is transversal to 
all the tangents of edges of 7 at v. We will see later that the end result of the 
regularized operator is independent of that additional edge, even better, that 
functions on graphs with only divalent vertices are annihilated by H E (N). 
With this preparation we will assume from now on that all vertices are at 
least trivalent. Numerate all the edges ej of the (so possibly extended) graph 
by some index I,J,K,.... Also we take all edges to be outgoing at a vertex 
as follows : by definition an analytical edge is an analytical embedding of a 
compact interval into E and a vertex v is a point of 7 such that there is no 
open neighbourhood U C E of v such that 7 D U is an embedded interval. 
By definition an edge is bounded by two vertices. Given an edge e of 7 with 
endpoints v, v' we subdivide it into two parts e = e o e' where e, e' respectively 
are outgoing at v, v' respectively. Note that e n e' is not a vertex of 7 so that 
e, e' are strictly speaking no edges any longer but we will continue to call them 
edges again as it simplifies the notation. 

1) Segments and arcs : 

Given an edge ej incident at a vertex v choose sj to be a segment of ej which 
is such that it 

a) is incident at v with outgoing orientation and 

b) which does not include the other endpoint vj of ej distinct from v. 

Also, given an unordered pair of edges ej, ej incident at a vertex v, let ajj be 
a curve which is such that that it 

i) intersects 7 in the endpoints of Sj, Sj distinct from v and 

ii) does not intersect 7 anywhere else. 

The diffeomorphism invariant properties of the position of the arc ajj will be 
specified more precisely below. We will see later that in order to meet the 
requirement that the constraint operator be symmetric on H we have to relax 
the prescriptions b),ii) above. 

2) Tetrahedra saturating a vertex : 

For each vertex v of 7 and each unordered triple of mutually distinct edges 
(e/, ej, ck) incident at v (the number of such triples is given by E{v) = 
n(v)(n(v) — l)(n(v) — 2)/6 where n(v) is the valence of the vertex) construct 
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eight tetrahedra saturating v as follows : 
We have a map 

(si, sj, s K ] au, a JK , a KI ) -> (si(A), s 2 (A), s 3 (A); ai 2 (A), 023(A), a 3 i(A)) where 
the labelling is such that the orientation of the the tangents at v is positive and 
we have indexed these six segments by the obvious tetrahedron A that they 
form. Choose the basepoint of A to be the vertex of 7 under consideration, 
that is, v (A) = v. 

We are now ready to construct the eight tetrahedra saturating v for the triple 

(e I ,e J ,e K ) : 

Let [0, 1] — > Sj(A)(t) and [0, 1] — > a i3 -(A)(i) be parameterizations of Sj(A) and 

Ojj(A) respectively. Define their "mirror images" by 

s,(A) a (t) :=2u - Si (A) (*), 

a Tj (A) a (t) :=2v a -a ij (A) a (t), 

a^(A)(t) := a^(A) a (t) - 2f[v - Sj (A)(l)]°, 

ay(A)(0 := ay(A)»(f) + 2t[ W - Sj (A)(l)] a 

respectively where w a are the coordinates of the vertex v. Here we have as- 
sumed that all objects lie in a chart containing v which is always possible 
by choosing the basic quantities Sj(A), £iy(A) small enough. Notice that by 
definition S;(A)(0) = v,S;(A)(l) = a ii (A)(0), s 3 - (A)(1) = a^(A)(l) so that 
S? (A)(0) = u,Si(A)(l) = a ?J (A)(0) = % -(A)(0), sj(A)(l) = a^(A)(l) = 
a i j(A)(l),s i (A)(l) = a fj ( A) (0),Sj (A)(1) = a^(A)(l). We did not make use 
of any background metric. We now form loops ajj, ay, aq and combine them 
with Sfc, to form seven more right oriented tetrahedra. Together with A 
these are the eight tetrahedra that we looked for. We will see that the seven 
"mirror" tetrahedra do not play any role at the end of the day so that the 
choice of adapted frame is irrelevant to define them. 

Although we will only be concerned later with the one-dimensional edges of 
the tetrahedra constructed, we will need also their two- and three-dimensional 
properties in an intermediate step : 

Choose any surfaces bounded by the edges of these eight tetrahedra and define 
the closed subset of X bounded by those faces of a tetrahedron A to be the 
region assigned to A. The only property of the region assigned to A which 
will be important is that A and its mirror images saturate v which is a diffeo- 
morphism invariant property and therefore everything will be independent of 
these regions. 

3) Diffeomorphism invariant prescription for the position of the arcs Ojj(A) : 
The following lemma, whose proof to the best of our knowledge is unpublished, 
shows that one can always choose two curves to lie in the x/y plane of an 
adapted coordinate system in which they take a standard form. 

Lemma 3.1 Let s 1 ,s 2 be two distinct analytic curves which intersect only in 
their starting point v. There exist parameterizations of these curves, a number 
5 > and an analytic diffeomorphism such that in the corresponding frame 
the curves are given by 

a) Si(t) = (t, 0,0), s 2 (t) = (0, t,0), t G [0,5] if their tangents are linearly 
independent at v and 
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b) Si(t) = (£, 0, 0), s 2 (£) = (£, £ n ,0), t G [0,5] for some n>2 if their tangents 
are co-linear at v. 

We will call the associated frame a frame adapted to si,s 2 . 
Proof : 

Given a frame, denote by 6j, i = 1, 2, 3 the standard vector with entry 1 at the 
i-th index and zero otherwise. 

First we show that any curve s can be mapped into a straight line by an 
analytic diffeomorphism. To that end, let us expand s(t) = f l {t)bi where 
f % are analytic functions of t. Since s is nowhere vanishing, at least one 
of the functions, say f 1 , has the property / 1 (0) 7^ and so it does not in 
a neighbourhood of 0. Choose b[ := — 62^1 _ ^3,^1 — 62 — 63 and 

{P\ Z 3 ') = (P, J 3 ), (P + f\ P), (P, P + P), (P + f\ / 3 + P) whenever 
(/ 2 (0), / 3 (0)) are 0, ^ 0), (= 0, ± 0), 0, = 0), (= 0, = 0). We conclude 
that we can write sit) = P(t)b\ + (P)'(t)b 2 + (P)'(t)b 3 =: g^t)^ where b\ 
form a basis and c/i(0) 7^ 0. It follows by the inverse function theorem that 
the equation x l = g l (t) can be inverted in a neighbourhood of and that 
{g l )~ l {x l ) is analytic in a neighbourhood of because g l (t) is of order o(t). 
We now see that the following diffeomorphism x %t (x 1 , x 2 , x 3 ) := is 
analytic and maps s(t) to s'(t) = x'(s(t)) = t(b[ + b' 2 + 63). Upon performing 
a constant diffeomorphism (that is, a GL{3) transformation) we can achieve 
that s'{t) = tbi. 

So we can assume that we have already mapped si, s 2 so that si(i) = £61. Now 
consider s 2 (t) = /*(£)6j. 

Case a) Since Si(0), £2(0) are not co-linear (which is a diffeomorphism invari- 
ant statement) it follows that not both of / 2 (0),/ 3 (0) can vanish. So let us 
assume that for instance / 2 (0) 7^ 0. By a similar argument as above we can 
write S2(£) = g l it)b\ + g 2 (t)b' 2 + g 3 {t)b 3 where 61, b' 2 , 63 are linearly independent 
and g*(0) 7^ 0. The analytic diffeomorphism x 1 ' (x 1 , x 2 , x 3 ) = (g l )~ 1 (x' t ) maps 
s'i(£) = x'(s(t)) = g v (t)bi, s' 2 (t) = t(bi +62 + 63). Now, a change of parameter- 
ization t' = g v (t) for Si and a final GL(3) transformation proves the assertion. 
Case b) Since Si(0),s 2 (0) are co-linear it follows that / 1 (0) 7^ 0,/ 2 (0) = 
/ 3 (0) = 0. Let n > 2 be the smaller of the two numbers n 2 ,ns defined by 
Pit) = o(£ n2 ),/ 3 (£) = o(t n3 ) (n is finite because the curves are not identical). 
Without loss of generality we may assume n 2 = n and now it follows by an 
already familiar argument that we can write s 2 {t) = g l it)bi + ^ 2 (£)6 2 + g 3 (t)6 3 
where 61, 6' 2 , 6 3 are linearly independent and (^(O), ig 2 Y n \0), (g 3 Y n \0) 7^ 0. It 
follows that g 2 (£) =: ^ 2 (s),^ 3 (£) =: ^ 3 (s) are analytic functions of the ana- 
lytic coordinate s := t n and that they are invertible in a neighbourhood of 0. 
We now define the analytic diffeomorphism 3 ) := ig l ) 1 ix l ) where 

g 1 = g 1 which maps s[(t) = 6 1 (g 1 )~ 1 (£), s' 2 (t) = bit + £ n (6' 2 + 63) and again a 
reparameterization for Si and a GL(3) transformation proves the assertion. 
□ 

Using this lemma we will now describe more precisely the choice of the arcs 
dij(A). We follow and extend (to the case of a pair of edges with co-linear tan- 
gents at their intersection) an elegant prescription due to Jurek Lewandowski 
|p~5f to which one is driven quite naturally. For the sake of self-containedness 
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of the present paper we repeat the argument here to the extent we need it. 
The discussion is rather technical and lengthy and the reader not interested 
in the details may skip the rest of the present item and just should assume 
that there exists a diffeomorphism invariant prescription of the topology of the 
routing. 

Let si, S2 be two segments of edges e±, e 2 of 7, incident at the vertex v. Their 
other endpoints t>i,i>2 are connected by an arc a. Basically we wish to avoid 
that the arc a intersects any other point of 7 different from i>i,t>2. Clearly, by 
choosing a to lie in a small enough neighbourhood of v the only danger is that 
a can intersect some other edge e different from e±, e 2 , also incident at v. With 
rising valence of v the number of topologically different possibilities of routing 
a between the edges e incident at v becomes rather complex so that we need 
to give a diffeomorphism invariant description of the choice of routing. 
By choosing a frame adapted to s\, s 2 we have partially fixed the gauge freedom 
associated with analyticity preserving diffeomorphisms. Consider any other 
frame (x lf ) connected to the identity (so that in particular x l ' xi > 0) and 
preserving the adaptedness of the frame. Thus x l '(s 1 (t)) = (t, 0,0) implies 
that, upon taking arbitrary derivatives with respect to t, that d™y' = d™z' = 
for all m > at v. Likewise, x l> (si(t)) = (0,t, 0) implies &^x' = d™z' = 
for all m > at v. However, for x l '(t,t n ,0) = (t,t n ,0) it is more subtle to 
characterize the partial derivatives. Consider the k— th derivative at t — of 
the function z'(t, t n , 0) = 0. Let k = mn+r where r = 0, .., n—\. It follows that 
we get a linear combination, with positive coefficients, of the following partial 
derivatives at v : d r x +ln d™~ l z' for I — 0, .., m which has to vanish. In particular, 
the case m = 1, r = 0, .., n — 1 shows that we get a linear combination of d r x +n z' 
and d r x d y z' so that d r x d y z' — at v for all r = 0, .., n — 1. The case m = 2, r = 
gives a linear combination of d x n z' = 0, d x d y z', d 2 y z' and the latter two terms 
do not need to vanish separately. Working out the precise coefficients of that 
linear combination we find that 



It turns out that this is all we need to prescribe the routing of the arc. 
The idea is that we wish to choose the arc a to lie in the coordinate plane 
bounded by si,S2 (that is, either in the quadrant x, y > or the wedge 
{(x, y,z); < y < x n , z = 0} respectively) and the question is under which 
circumstances we can guarantee that then a does not intersect e. 

The easiest case is that ^(0) 7^ 0. We show that sgn(e 2; (0)) is preserved under 
change of adapted frames. Namely e z '(0) = z' a e a (0) = z' z e z (0). By analyticity 
it then follows that e is curved away from the x/y plane in a neighbourhood 
of v and so does not intersect a for no choice of adapted frame and therefore 
for no frame. 

The case e z (0) = is more tricky. Without loss of generality we may assume 
that at least e^O) 7^ (otherwise e(0) = 0; switch the role of x and y if 
necessary for the case that si, s 2 do not have co-linear tangents at v. If they 
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do have co-linear tangents at v and e^O) = then e y (0) 7^ and it follows 
that a segment incident at v of the projection of the curve e into the x/y plane 
lies above the parabola y = x n which describes s 2 and so a cannot intersect e 
anyway) . 

We can make a distinction between two situations. 

Situation A : There exists a finite number m' such that (e*" 1 )) 2 (0) 7^ and m > 
m', m < 00 such that (e^) y (0) ^ are the first non-vanishing derivatives. 
The combination m' = m = 00 is excluded as otherwise Si and e would overlap 
in a finite segment (here we have used the analyticity of the edges). We then 
readily verify that under a change of adapted frame we have that the first 
non-vanishing derivative of the z-component of e is given by (e*™)) 2 (0) = 
z' tZ (e^ m '^) z (0) so that the sign of (e( m ')) 2 (0) is again preserved and e is curved 
away from the x/y plane so as not to intersect a. 

Situation B : There exists a finite number m such that (e (m) ) ? '(0) ^ and 
m' > m, m' < 00 such that (e*" 1 )) 2 7^ are the first non- vanishing derivatives. 
Again m = m' = 00 is an excluded possibility. 

Case that Si(0),S2(0) are not co-linear. 

If e does not point into an octant of the frame where both x/y are positive then 
again there is no danger that a intersects e since we choose a to lie in the x/y 
plane with positive x/y components as said above. Since x' x ,y' y 7^ 0,d x y' = 
d k y x' = OV/c it follows that the sign of the first non-vanishing derivatives of 
e x , e y are preserved under changes of the adapted frame. 
So let us assume that e does point into an octant where both x/y are positive. 
This means that e x (0), (e( m ))^(0) > 0. 

The first possibly non-vanishing derivative of the ^-component of e is given 
by (e( m+1 )) 2 '(0) = z' xy e x (0)(e^) y (0) + ^ 2 (e( m+1 )) 2 (0) and we see that this can 
have any sign for a suitable choice of z'. We use this freedom to further fix 
the frame such that this sign is positive. Choose z'(z,y,z) = z + f3xy. This 
satisfies all requirements on z' at x = y = z = and we see that upon choosing 
one and the same f3(s\, s 2 ) large enough we can beat the terms z' 2 (e (m ~ fl )) 2 (0) 
for an arbitrary (but finite) number of edges e as to make (e(™ +1 )) 2 '(0) > 0. 

Case that Si(0),s 2 (0) are colinear. 

If e does not point into the wedge y < x n ,x > then again the arc a cannot 
possibly intersect e if we choose it close enough to the vertex. This time we 
only have that x' x , y' y 7^ 0, d x y'Wk which implies only that the sign of the first 
non- vanishing derivative of e y is preserved, we even find that x' x — y f , y — 1 
so that even its value is preserved. Now, if m > 2 then the sign of e x (0) is 
preserved. If m = 1 then it is not necessarily preserved but then it is true 
that, since n > 2, the projection of e into the x/y plane, which lay outside 
the wedge in a small enough neighbourhood of the vertex, is still outside the 
wedge. Therefore the condition that e does or does not point into the wedge 
is preserved. So let us assume that e points into the wedge, in particular, 
e x (0), (e( m ^(0) > 0. 

We notice that the first derivative of e 2 at t — which involves A := d x d y z' 
at v is of order n + m and given by a term CA^4(e x (0)) ra (e^) y (0), ca > 0. 
Likewise, the first derivative which involves B := d y z' at v is of order 2m and 
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given by a term c B B[(e^) y ($)f , c B > 0. 
Subcase I) w! < min(2m, m + n) 

In this case the first non-vanishing derivative of e z at t — is of order ml and 
involves a term z' z (e^ mn> ) z (0) . By choosing the coefficient z' z large and positive 
enough we can beat any possible contribution involving A, B and preserve the 
sign of (e (m,) ) 2 (0). 
Subcase II) to/ > min(2m, m + n) 

The first non- vanishing derivative of e 2 at t — is of order m + n if n < m and 
involves only the term proportional to A, of order 2m if n > m and involves 
only a term proportional to B if n > m and for n = m the derivative is 
precisely given by 

( e ( 2n )) 2 '(0) = B(e {n) ) y (0) ^ 2n ~ 1 ^ [(e w ) 1 '(0)-n!(e 1! (0)) fl ] 

where we have made use of the relation between A, B as stated above. By 
choosing e x (t) as a parameter we may assume that e x (t) = t. Then we may as- 
sume (e^) v (0) < n\ : Namely, if (e^) y (0) > n\ then again a segment incident 
at v of the projection of the curve e into the x/y plane, whose y— component 
looks like e y (t) — t n + o(t n+1 ), lies above the parabola y = x n which de- 
scribes s 2 and cannot possibly intersect the arc a. If (^"^(O) = n\ then wc 
must have that either e y has another higher order non-vanishing derivative 
(e(™)) v (0) 7^ 0, to > n or, if that is not the case, m! < oo for otherwise e would 
coincide with s 2 which we excluded. In the latter case, the first non-vanishing 
derivative of e z at t — would be again z' z (e^ m ~ ) ) z (0) because by definition 
of a diffeomorphism preserving the adaptedness of the frame, all derivatives 
of z'(e(t)) = z'(t,t n ,e z (t)) at t — which do not involve at least one partial 
derivative with respect to z must vanish. So again the sign of (e( m '^) 2 (0) would 
be preserved. In the former case, the first non-vanishing derivative of e z at 
t — and which is proportional to d^.d l y for some k > and some I > at v 
is of order m + n and proportional to d"dyz'(e^) y (0) which follows from the 
fact that all contributions which do not contain at least one factor of (e( m ))^(0) 
must vanish due to the adaptedness of the frame. Now we are back to either 
m' < n + m = min(m + n, 2m) or m' > m + n and we have the case m > n. 
It follows then that if we choose any positive number f3 and some large enough 
positive number a > and the diffeomorphism (which satisfies the condition 
between A, B as given above) z' := az + j3(y 2 /2 — 2yx n /(n\) 2 ) then the im- 
age of the edge e under this diffeomorphism will be such that it is curved 
away from the x/y plane as before for m' < min(m + n, 2m) (where it is 
understood that we take m to be the next to leading order of e y (t), with 
positive coeeficient, in case that e x (t),e y (t) — t n + o(t m ) and m = oo other- 
wise) and for m' > min(m + n, 2m) it is curved into the upper half space if 
n > m, into the lower half space if n < m and into the lower half space if 
e x (t) = t, e y (t) = kt n + o(t n+1 ), k < 1 (for the other cases there is no routing 
to be chosen). This can be achieved for an arbitrary number of edges e with 
the same a(si, S2), f3(s±, S2)- Notice that the topology of the routing is dif- 
feomorphism invariantly defined as the numbers n, to, to' are diffeomorphism 
invariant. 
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Concluding, we choose an adapted frame and a small enough neighbourhood 
of v and an arc ciy(A) going through that neighbourhood, which lies in the 
plane bounded by s«(A), Sj(A) such that the routing through all other edges 
incident at v is the one described above. 



4) Tetrahedra away from the vertices : 

Denote by D(A) the closed region in E filled by the eight tetrahedra con- 
structed from a triple ej,ej,ex as outlined in 1). Also consider their union 
D(v) = U u (a)=«-D(A) and their complement with respect to D(v), that is, 
D(A) := D(v) — D(A). We triangulate D(A) as outlined in 1) by the eight 
tetrahedra constructed. The sets D := S — \J v ^y^D{v) and D(A) are trian- 
gulated arbitrarily. As we have argued above, the final result will not depend 
on these tetrahedra because they do not intersect a vertex of 7. This follows 
from the fact that all the tetrahedra different from the ones constructed in 1) 
have a basepoint different from any of the vertices of 7 since by construction 
the tetrahedra described in 1) saturate them. 

5) Closeness to a triangle : 

Having fixed the topology of the routing of the arcs as in 2) we now choose, in 
the standard frame, the arc fly (A) to be as straight as possible so that ot^ (A) 
looks like a triangle, as much as the routing allows. This will then justify the 
approximation h a ..( A) « 1 + ^s i (A) a (0)s j (A) 6 (0)F ab (i;(A)) where [0,St] is 
the subinterval of [0, 1] corresponding to s»(A) as compared to the whole edge 
of 7 of which it is a segment. 



3.2 Final Regularization of the Euclidean constraint 



Let us now write the integral over S (SO) for the classical theory as follows 

x-u vev(l) D(v) + Jd(v) ~ Id + ug ^ 7) E(v) v ^l}L(a) + L(J 

-^L+^^^l^ / A ,+ E Jj- (3-6) 

In the classical expression (|3.3|) we can replace, for each A, ^A'sDfA) Ja' by 8 / A as 
all the tetrahedra shrink to their basepoints so that the dependence on the mirror 
images of A drops out as promised. 

Now we quantize the resulting expression ( |3.3|) as outlined above and it follows from 
the considerations in section 3.1.1 that on the cylindrical subspace labelled by the 
graph 7 the action of the Euclidean Hamiltonian constraint reduces to 

H% 7) [N\f := Hf[N\f = £ N v ^- £ Hif=: £ N v H*f (3.7) 

veV{-y) -^V v ) v(A)=v ueV(7) 

This furnishes the regularization step. In particular, the expression ( |3.7| ) is finite (i.e. 
a cylindrical function) independently of T because the number E(v) is determined 
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by the graph and not by the triangulation and so does not change as we make the 
triangulation finer and finer. 

Let us now show that divalent vertices do not contribute : if v is a divalent vertex 
of 7, if S3 is the segment of the edge added in requirement 0) and if au is the 
corresponding loop along segments s±, S2 of the two edges ei,e2 of 7 incident at v 
then 

16 

H* = -^tr((V 2 - Kl 2 )h s ,[K^ V}) . (3.8) 

The terms involving a 2 3 , 0:31 drop out which follows from the fact that the volume 
operator annihilates divalent vertices as is obvious from the expression ( |2.8| ). Notice 
that if s 3 was not transversal to s 1; s 2 at v then (|3.8| ) would vanish trivially because 
the volume operator annihilates vertices which are such that all incident edges have 
co-planar tangents. But even so, ( |3.8|) vanishes : since the end result of applying 
( |3.8| ) to a gauge invariant function / cylindrical with respect to a graph 7 must be 
gauge invariant, it cannot depend on S3 (to see this, expand H E f into spin-network 
states. Each of these spin-network states can colour S3 only with spin because 
there is no edge of 7 available in order to combine with the other endpoint of S3 
different from v in a gauge invariant way). Moreover, it is easy to see that only a 
term proportional to h S3 Vh~ 3 l survives. Now from the fact that the space of vertex 
contractors for divalent vertices is one- dimensional we see that h S3 Vh~ 3 f = \f so 
that tr([h ai2 — h ai2 ]h S3 Vh~ 3 )f = Xti([h ai2 — h ai2 ])f = where we used the SU{2) 
Mandelstam identity. The argument actually extends to the case of vertices of ar- 
bitrary valence but such that all incident edges have co-planar tangents. So we see 
already that functions cylindrical on graphs with only such degenerate vertices are 
annihilated by H E (N). This is a feature which is shared with previous regulariza- 



tions 28, m 



It is amazing that one got expression ( |3.7| ) almost for free once one knows that 
the volume operator is well-defined on holonomies, no ill-defined products of dis- 
tributions arise, we do not encounter any singularities, no renormalization of the 
operator is necessary. Note that we have no problems in ordering V to the left or 
to the right of the holonomies involved as V has a finite action on holonomies of A 
as is clear from 



3.3 Cylindrical consistency 

We have now produced an uncountable family of Euclidean Hamiltonian constraint 
operators (H E (N))^, one operator for each graph 7 given in ( |3.7[) . What we need 
to make sure is that these operators are the projections to cylindrical subspaces of 
one and the same operator on 7i. This requirement will lead to some modifications 
of ( |3.7| ) (while keeping the classical limit to be still H E (N)). 

Also we would like to construct one version of H E (N) which is symmetric. The way 
it stands, not even the projections of H E (N) given in ( |3.7| ) are symmetric operators 
on 7i. Therefore, for the symmetric operator, we first order each term associated 
with a tetrahedron symmetrically . The result is (using that V is symmetric on 
H and h\ = ~h e = {h^)' 1 on H) that, in case we wish to construct a symmetric 
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operator, we replace ( |3.7| ) by 



h*[N]f:=j:NXf,h*=— £ h%, 

v£l ^\ V ) v(A)=v 



h 1 , 



3^^Hr({^. (A) , h Sh[A) [h; k \ A) , V]}) (3.9) 



where {., .} denotes the anti-commutator, while we stick with ( p.7|) if we do not wish 
to construct a symmetric operator. With this choice, each h®[N] in ( |3.9|) separately 
becomes a symmetric operator on 7i. Notice that this is far from guaranteeing that 
H E (N) itself is symmetric (see |1]]). Both operators in ( |3.7| ), ( |3.9| ) have the dense 



domain Cyl (A/G) inherited from V [26|, the thrice differentiable cylindrical func- 



tions on A/G- 

We now come to make both (|3.7|) , (|3~9| ) cylindrically consistent. Let A e = X l {h e )X l (h e 
be the Casimir operator associated with an edge e of 7 incident at v. Note that 
since X(h e ) is right-invariant we have X(h e ) = X(h s ) for any segment s of e as 
long as s is incident at the outgoing endpoint of e. We can now define an edge 



projector p e := 8(j e ) where j e := Jl/4 — A e — 1/2 has spectrum 0,1/2,1,3/2,.. 
and 9 is some smooth function on 11 which vanishes on (—00, 1/8]) and equals 1 on 
[3/8,oo). The operators p e are all commuting among each other and symmetric. 
The effect of this operator when applied to a function / cylindrical with respect to a 
graph 7 is to annihilate that function if 7 and e do not intersect in a finite segment 
incident at the outgoing endpoint of e and to leave it invariant otherwise. From 
these projectors we construct a tetrahedron projector p& '■= p Sl (A)Ps 2 (A)Ps 3 (A) and 
a vertex operator E{y) := J2v(A)=vPa- We then define the following self-consistent 
(up to a diffeomorphism) family of (symmetric) operators 



^Jfe for (3.7) 

PA j^E PA 

ueV(-y) v(A)=v y y/E(v) A y/E(v) 



H*[N):= £ 8N V J2 { **\ pa ' for q (3-10) 



We note that ( p,10|) still has the correct classical limit : up to terms of higher order 
in h its action on cylindrical functions is still given by ( J3.71 ) which was shown to 
have the correct classical limit. 

The symmetry of each member of the family defined in the second line of ( |3.10| ) is 



obvious. The self-consistency of both families defined in ( |3.10| ) can be checked as 
follows : 

A graph 7 C 7' can be obtained from 7' by a finite sequence of steps consisting of the 
following two basic ones : a) delete an edge of 7' and b) after removing an edge e, if 
one (or both) of the former endpoints v of e is now divalent and v — ex l~l e?, where 
e±, 62 are analytic extensions of each other then combine e±2 ■= e\ o 1 to an edge of 
7, that is, delete a vertex. In case a) it is clear that each term in ( |3.10| ) corresponding 
to a tetrahedron which involves the removed edge vanishes identically while in case 
b) all terms corresponding to the removed vertex vanish when applied to a function / 
cylindrical with respect to the graph 7. Moreover, E(v ) reduces to the correct value 
on /. Finally, it follows from our manifestly diffeomorphism-invariant prescription 
of a loop-assignment that Hy f and H^f are related by a diffeomorphism. In more 
detail, we have the following : if 7' is bigger than 7 and if e[, e' 2 are edges of 7' which 
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are the parts of the edges ei,e 2 of 7 incident at the same vertex v at which ei,e 2 
are incident then the loop made from the corresponding segments s^s^ and S\, s 2 
are diffeomorphic thus guaranteeing cylindrical consistency up to a diffeomorphism. 
This is enough to show consistency. 

In the sequel we will prove that both operators in (|3.10|) share the properties of 
diffeomorphism covariance and anomaly-freeness. To see that the second operator 
in ( |3.10| ) is actually symmetric requires a modification of the regularization which 
does not spoil those properties. We will come back to the modification in [|TJ. 



3.4 Diffeomorphism covariance 

According to the programme of algebraic quantization proposed in [[13] the solutions 
to the Euclidean Hamiltonian constraint are diffeomorphism invariant distributions 
if) on $ := Cj\°°(A/Q) such that for all lapse functions N 

ty[H E [N](f)] = V G $ . (3.11) 

Now take = / to be any function, cylindrical with respect to some graph 7, in the 
domain of H E then (|3.11| ) amounts to J2vev(~/) ^[Hyf] = f° r an N v ,v G V(j) and 



therefore we find that we need to satisfy 



9[H?f>r) = V 7, A G Cjt(A/G), v G V( 7 ) ■ (3.12) 



Equation ( |3.12j ) is actually quite unexpected : in section 3.1.3 we formulated the 



requirement of diffeomorphism-covariance in terms of the constraint at a point and 
is was far from clear that such an operator actually makes sense. Equation ( p. 12 ) is 
the precise formulation of that concept and it is manifestly well-defined. 
Our triangulation adapted to a graph was geared at being diffeomorphism covari- 
ant for each of its vertices separately meaning that each of the operators H E and 
even each of the operators p&, E(v), H E , separately is covariantly defined. For 
p&,E(v) this follows from the manifest covariance of p e , namely U(ip)p e ll(tf)~ 1 = 
p v ( e ) for any if G Diff(S) which in turn is a consequence of the fact that p e is 
defined in terms of h e (A). For the operators H^yh^ we argue as follows : let f 
be cylindrical with respect to a graph 7 and let if G Diff(E). Then the tetrahe- 
dra A((f(j)), ip(A(j)) are not necessarily equal to each other. However, the graphs 
7,7' = ^(7) are diffeomorphic, therefore the topology of the routing of the arcs 
ay(A(7)) through the edges of 7 and of the arcs Oy(A(7')) through the edges of 7' 
as specified in section 3.1.3 is the same since that prescription depended only on the 
topology of the graph. More specifically, this prescription was shown to be indepen- 
dent of the frame and therefore coincides for any two vertices v , v' of graphs 7, 7' 
for which neighbourhoods U, U' exist such that U D 7 and U' H 7' are diffeomorphic. 
Now just choose a diffeomorphism f' such that ^'(7') = 7' and such that y? / (A(7 / )) = 
if(A('~f)) for that specific A(7) (notice that ( |3. 12|) is a linear combination of terms, 



each of which depends on only one specific A so that we can adapt if' to A). Such 
a diffeomorphism clearly exists : there are diffeomorphisms which leave the image 
of the graph invariant while moving its points and off the graph it can put the 
arc dy(A) into any diffeomorphic shape. It follows from these considerations that 
V \y')h^i , y ylJ (f 1 )' 1 = ^(A( 7 ) which is what we wanted to show. 
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Then, obviously, the number ty[H E f y ] depends only on the diffeomorphism class 
of the loop assignments ofy(A) (this was first observed in fll4l). Therefore, in this 
diffeomorphism invariant context, the loops ctij(A) can be chosen as "small" and 
the triangulation as "fine" as we wish, the value of (|3.12| ) remains invariant and in 



that sense the continuum limit has already been taken. 

Diffeomorphism covariance is therefore a sufficient requirement for our quantum the- 
ory to correspond to a continuum theory. 

One might wonder what happens if one actually takes the limit and sends A — > v (A). 
It is easy to see that the result vanishes trivially which is not what we want. This 
happens due to the fact that after applying H E (N) to a cylindrical function only a 
finite number of terms survive : since a cylindrical function is already determined 
on smooth connections we can, in the limit, actually replace Hj[ by an integral over 
A as in (|3.1|) but the limit corresponds to a point which has zero Lebesgue measure. 
The fact that the limit A — ► v(A) is trivial is strange at first sight because one is 
used, from the lattice regularization of, say, A</> 4 theory, that the continuum the- 
ory is only recovered if we take the lattice spacing to zero, that is, one takes a 
continuous cut-off parameter to its continuum value. In our regularization such a 
parameter simply does not exist and the reason for that is the underlying diffeo- 
morphism invariance of the theory. This shows that the limit A — > v(A) is in fact 
inappropriate. 

3.5 Anomaly-freeness 

Although our assignment is covariant and therefore the continuum limit is already 
taken (in the sense explained above) so that it seems that the regulator is entirely 
removed, the operator (|3.10|) still carries a sign of the regularization procedure : 
it depends on the diffeomorphism class [T] of the triangulation assignment which 
labels the freedom that we have in our regularization scheme. It is therefore not an 
entirely trivial task to check whether our operator H E [N] is anomaly-free, meaning 
that [H E [M], H E [N]]f vanishes for any cylindrical function / and lapse functions 
M, N when evaluated on a diffeomorphism invariant state tp. The reason why we do 
not check the commutator on tp immediately is because ip is a distribution jn| and 



so does not lie in the domain of (H (M)y . Therefore, the formal anomaly-freeness 
condition [(H E (N)y , (H E (M)y]ip = has to be interpreted in the usual weak sense 
([(H E (N)Y, (H E (M)Y]^)(f) = iP([H E (M),H E (N)]f) = for each test function /, 
that is, every smooth cylindrical function / £ $. 

If the theory is anomaly-free, then, since ip is a generalized eigenvector |13|] for the 



exponentiated diffeomorphism constraint {/(</>), £ Diff(S) with eigenvalue 1, we 
expect that in the last equation the argument of if} is identically zero if some finite 
diffeomorphisms can be removed. This expectation turns out to be precisely correct. 
We show now that a solution to the anomaly freedom problem is obtained 

a) for the non-symmetric operator in case we attach the edges ay(A) irrespective of 
their differentiability at their endpoints with respect to s^(A),s 3 -(A) as prescribed 
in section 3.1.3 

b) for the symmetric operator only if all the loops otij(A) are chosen to be kinks with 
vertex at v(A) ! That is, the arc ay (A) joins the endpoints of ej(A), Cj(A) in at least 
a C 1 fashion (see ffl for the details of the attachment). An arbitrary attachment of 
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ay (A) is insufficient to guarantee anomaly- freeness. 

Consider for simplicity a graph 7 which only has one vertex v and that it is two- 
valent, for instance a kink (we ignore for the moment that functions on such graphs 
are actually annihilated in order not to veil the argument, the problem shows up 
on higher valent graphs). Acting once with the Hamiltonian constraint on a func- 
tion / cylindrical with respect to 7 we get a function cylindrical with respect to 
a graph 7' which contains 7 and an additional edge e which intersects 7 in ver- 
tices Vi 1 v 2 and such that the tangents of e, 7 at these new vertices are, a priori, 
linearly independent. Acting once again with the Hamiltonian constraint it acts 
now non-trivially at all three vertices. Therefore, the commutator becomes now 
schematically [H E (M), H E (N)]f = (M V H E + M Vl H E + M V2 H E )N V H E f - (N V H E + 
N Vl H E + N V2 H E )M v H E f = (M, 1 ^-^ 1 M„)^^f/ + (M, 2 ^-^ 2 M,)^Ff/ 
which does not manifestly vanish even if [H E , H E ] = for v 7^ v '. One sees that 
what has to be avoided is that the Hamiltonian has non-trivial action at Vi, v-i. The 
idea is to exploit that the volume operator acts trivially on vertices which are such 
that the tangents of all edges incident at it lie in a common plane. 

It is here where the volume operator as defined in the second reference of is 
selected by the dynamics of the theory while the operator as defined in /]|3J/ has to be 
rejected (if one follows the approach advertised here) ! 



Namely, the volume operator p6[ does not annihilate co-planar vertices. It is 
very appealing that the requirement of anomaly-freeness provides us with a selection 
rule among the operators |[25|| , p6| which, a priori, from a purely kinematical point 
of view, are both bona fide quantizations of the classical volume functional. 

Now, for the non-symmetric operator there is no problem at all because by 
construction of the triangulation assignment the endpoints of ay (A) form always 
new three- valent vertices of 7 U ay (A) but the edges incident at them only have 
two independent tangent directions there, namely those of e$(A) and ay (A). But 
for the symmetric operator we need to adjoin the smooth exceptional edges (see |L]]) 
in at least a C 1 fashion because all the segments ay (A) which intersect the same 
edge e of the skeleton of the given graph at all, intersect it at the same point and 
in general for higher than valence two the volume operator is non-vanishing on (not 
necessarily gauge-invariant) cylindrical functions unless the tangents of triples of 
edges at vertices are linearly dependent. 



In [14, 15] the loop assigned does not have the topology of a kink but, the topology of 



a triangle. But because these operators do not involve the volume operator, rather 
they depend on an operator corresponding to e a b c €ijkE° 'Ej which does not vanish if 
there are only two independent tangent directions in the problem, we expect these 
operators to be anomalous. However, if the loop assigned would have the topology 
of a kink as well we believe that the anomaly could be removed once these operators 
are rigorously defined. 

Theorem 3.1 The Euclidean Hamiltonian operator H E (N) = (H E (N), D 7 ) as de- 
fined by \3.1(\) is non- anomalous. 



Proof : 

Let / be a function cylindrical with respect to a graph 7 and let A (7) denote 
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the various tetrahedra attached to it in applying H E (N). We may, without loss 
of generality, assume that / is a spin-network function. Clearly the functions 
^A( 7 )/j ^A(j)f depend on the graph 7 U A(7). We are being very explicit here 
in the dependence of the tetrahedra on the graph because this will be essential 
in what follows. We may assume that all the tetrahedron projectors p& are non- 
vanishing on /, that is, the dependence of / on all the edges of 7 is non-trivial. 
Then E(v)f = E{v 1 r ))f = n(n — l)(n — 2)/6/ where n is the valence of v in the 
graph 7. Now consider A' with v' := v(A') 7^ u(A) =: v. Then pA'H E f = H E f 
since the segments of edges of 7 which may have been removed in H^f are at a 
vertex different from the vertex v(A'). It follows that E(v')H E f = E(v',j)H E f. 
If, however, v — v' then it is possible that pA'9 = where g is a spin-network state 
appearing in the decomposition of Hj[f. Now, for the symmetric operator, since 
we project with p& before and after applying h&, such g are automatically removed 
so that either g = or E(v')g = E(v,^f)g again. For the non-symmetric operator 
we get instead E(y')g = E(v,j,g)g. With this preparation we compute for the 
non-symmetric operator (it is understood that only terms with Pa( 7 uA(7))-^a(7) ^ ^ 
are kept) 



v<=V(-y) ^\ V iU v{A(-y))=v 
ue y( 7 ) V ' '/ v(A(j))=v «'eV(7UA(7) u(A( 7 uA(7)))=u' 



H E [M]H E [N]f = £ £ H E [M]Hl {l) f 

E ^T~\ £ M V #A( 7 UA( 7 )) £^#A( 7 )/ 

(3.13) 



and similar for the symmetric operator, where the essential step has been the last 
one where all the contributions from V (7 U A (7)) — 1^(7), where v{A{ r y)) G ^(7) have 
been removed due to the fact that all these vertices are co-planar (for the symmetric 
operator this holds because of the kink property of those vertices as explained in 
Lemma 3.1 and the subsequent Collorary). We now compute (|3.13 ) again with the 
roles of M, N interchanged, subtract from ( |3.13| ) and obtain 



[H E [M],H E [N]]f 

1 

v,v'eV(-y) V ' '/ t;(A(7))=v,«(A(7UA(7)))=v' 

E 



M V ,N V - M V N V , 



, . ^ , E(v)E(v') 

X E [-^A'( 7 UA( 7 ))-^A(7) _ ^A( 7 UA'(7))-^A'( 7 )]/ 

«(A(7))=v(A(7UA' (-y)))=v,v(A' ( 7 ))=v(A' (7UA(7)))=u' 

(3.14) 

where in the second step the notation v < v ' assumes that we have ordered the ver- 
tices of 7 somehow which allowed us to write a sum over unordered pairs of vertices. 
Also we used that the antisymmetric product of lapse functions vanishes at equal 
vertices which was crucial in replacing E(v') by E(v',j) =: E(v') so that we may 
imagine to absorb this number into the lapses, N v /E(v) —>■ N v (as explained above, 
this happens for the symmetric operator already before taking the commutator). 
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Formula ( 3.14J) is valid for the symmetric operator as well if we just replace Hj[ by 



hf. 

It is far from obvious whether fl3.14| ) vanishes or not. Indeed, for genuine 7 and gen- 
uine choices of loop assignments depending on the graph, ( |3.14|) is a non-vanishing 
cylindrical function of positive L2 norm. We now evaluate a diffeomorphism invari- 
ant state on (|3.14|) . We can take the antisymmetric product of the lapse functions 



evaluated at different vertices out of the integral over A/Q and it will be sufficient 
to show that 

V ; [(-^A'(7UA(7))-^A(7) — -^A( 7 UA'(7))-^A'( 7 ))/] = (3.15) 

for each choice of v(A(j)) = u(A(7 U A'(7))) = v,v(A'("/)) = v(A'(-y U A(7))) = v' 
separately. To see this, notice first that the members of the first pair of tetrahedra 
given by (A (7), A(7U A' (7))) as well as the member of the second pair of tetrahedra 
given by (A'(7), A'(7U A(7))) are diffeomorphic. This follows immediately from the 
fact that v ^ v' so that there are disjoint neighbourhoods U of v and U' of v' where U 
and U' respectively contain both members of the first and second pair of tetrahedra 
respectively. Let ip G Diff(S) be chosen such that 7 D U' is left invariant and but 
that A (7 U A'(7)) = <p(A(j)). Likewise, choose <p' G Diff(E) such that 7 n U is left 
invariant but that A'^UA^)) = <^'(A'(7)). Then, using diffeomorphism invariance 
of ip we find that the left hand side of the (|3.15|) becomes 

tp[(U(ip')H E , (y) H E h) - U(<p)H E (y) H E ^)f] 
= mi^Hi^Hl^U^- 1 ) - U^)Hl h) Hl l{l p{ V ')- l )f] 
= ##f, (7) ,tff (7) ]/]. (3.16) 

That is, we were able to "match" the tetrahedra using diffeomorphism invariance. In 
the second equality we used the invariance U(<p)f — f f — U^ 1 (ip)f and similar 
for ip' in order to write the commutator in such a way that it becomes manifestly 
antisymmetric if we replace Hj[ — > h^. 

Now we just need to use Collorary 3.1 of (the non-symmetric operator is treated 
in a comment after this collorary) to see that the commutator in ( |3.16| ) vanishes 
identically because the part V v of the volume operator involved in -£^(7) d° es n °t 
act on the holonomies along edges incident at v' involved in -£^'(7) an< ^ v ^ ce versa 
since the vertices v,v' are different. That completes the proof of anomaly-freeness. 
□ 

A couple of remarks are in order : 

• The classical constraint algebra is given by {H E [M], H E [N}} = / E (MiV \ a -N M ja )q ab V b , 
where V a is the classical diffeomorphism constraint and q ab is the inverse metric ten- 
sor. Naively, one would expect that the quantum version of that would be 

hfi E [M\,H E [N\] = J^(MN, a -NM, a )*q ab V b * (3.17) 

where the *...* is to indicate that the operators that appear have to be regularized 
and ordered appropriately. The immediate problem with (|3.17| ) as, widely discussed 



in the literature, is that the constraint algebra is not a proper algebra, the structure 
functions depend on the canonical variables through q ab and if in the quantum the- 
ory, following the Dirac approach, the generator of the diffeomorphism constraint 
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does not appear to the right of q ab in ( p. 17] ) then one would not expect a genuine 
diffeomorphism invariant state to be annihilated by the commutator of two Hamil- 
tonian constraints : One says that the constraint algebra has an anomaly and that 
the quantum theory does not correspond to the classical theory because any element 
of the kernel of both the diffeomorphism constraint and the Hamiltonian constraint 
would have to satisfy the additional requirement that the right hand side of (|3.17|) 
annihilates it which reduces the number of degrees of freedom more than the classi- 
cal theory would do. 

One can then ask the question whether there exists a consistent quantization (reg- 
ularization and factor ordering) of H E (N) such that there is no anomaly. In [^7| 
the authors investigate a wide class of finite-dimensional theories (gauge systems) 
with a Hamiltonian quadratic in the momenta and a constraint algebra which mim- 
ics (|3.171) algebraically in the sense that H is replaced by the generator of a gauge 
transformation and q ab is replaced by some non-constant function of the canonical 
configuration coordinates. The authors find that a consistent quantization can be 
obtained but never in such a way that the Hamiltonian constraint is a symmetric 
operator. The intuitive reason for this is clear : since the generator V a of the uni- 
tary representation of the Diffeomorphism group has to be self-adjoint and since the 
metric tensor should be at least symmetric as well, then the left hand side of ( |3.17|) 
should appear in the symmetric ordering q ab Vb + V a q ab (or an even more complicated 
symmetric operator) whenever H E (N) is a symmetric operator. Now, in one of its 
versions, we actually did order H E (N) symmetrically ! Why, then, is that not in 
contradiction to the arguments given in |27[ and references therein ? 
As is often the case with "no-go theorems" , one needs to carefully check the assump- 
tions. The assumptions underlying the considerations in |Z7|] , when applied to our 
case, contain at least the following list (we do not list the assumption that one only 
has a finite number of degrees of freedom since we do not want to blame the failure 
of the theorem on that) : 

1) The operators H E (N), q ab (x), V b ( x) can be regulated and densely defined on 7i. 

2) The classical Hamiltonian constraint is a bilinear form in the canonical momenta. 

3) The classical Diffeomorphism constraint is a linear form in the canonical mo- 
menta. 

Only assumption 3) is satisfied here, the rest is violated : 

1) As was shown in || the representation of (one-parameter subgroups of) the 
diffeomorphism group on 7i is not strongly continuous. Therefore, by Stone's the- 
orem, there does not exist a self-adjoint generator. Also a operator correspond- 
ing to q ab (x) is entirely meaningless in our representation : one could write it as 
q ab = E^E b / det(g), however, both nominator and denominator are meaningless 

0- 

2) The dependence of H E on the momenta is through {A a , V} which is not even an 
analytic function of Ef. 

We conclude that the considerations of |27| are not applicable in to our case. 
How then can we even define a theory to be anomaly-free if we cannot define the 
quantizations of the generators of the symmetries of the classical theory ? The an- 
swer is that the quantization of a classical theory does not force us to make sense 
out of every classical function as an operator. The most general question we can 
ask is the following : 
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a) Can we make sense out of the right hand side of ( |3.17| ), that is, can we make 
sense out of 

/ (MN a -NM, a )q- b V b ? (3.18) 

The answer is, trivially, in the affirmative : up to higher order in % the commutator 
[H E (M),H E (N)}/(ih)} reproduces J s (MN, a - NM, a )q ab V b so that we can just define 
( |3.18| ) by that commutator because by construction it is well-defined. 

b) The fact that the classical Poisson bracket vanishes on the constraint surface of 
the phase space defined by the diffeomorphism constraint translates in the quantum 
theory into the requirement that ( |3.18j ) should vanish on diffeomorphism invariant 
states and this we checked to be the case. 

To conclude, in our case the question of whether a factor ordering can be found such 
that V a stands to the right in ( |3.18| ) cannot even be asked. Therefore the quantum 
constraint algebra cannot have any close analogy with the classical constraint al- 
gebra. Given the fact that we can only define an operator corresponding to finite 
diffeomorphisms the structure of the constraint algebra that we do expect is pre- 
cisely the one displayed in ( |3.16|) , namely that the commutator on a state / is a 
sum of terms of the form U((p)ABf — U(<p')BAf where [A,B]f = and A, B are 
symmetric operators. The diffeomorphisms (p, <p' G Diff(E) are highly ambiguous 
since their action has only been specified on a finite graph induced from the one 
underlying /, in particular U(ip)f = U(ip')f = f. Using this ambiguity we can 
write the right hand side of the commutator [H E (M), H E (N)]f again in manifestly 
anti-symmetric form U (ip)ABfU ((p')^ 1 — U((p')BAU {ip)^ 1 f as in ( |3.16| ) so that there 
are no factor-ordering contradictions of the kind discussed in [27| at all ! 
• It will turn out that we do not need to apply any group averaging |Tj| with 
respect to the non-symmetric Hamiltonian constraint operator but rather can com- 
pute the solutions by direct methods. This is important because if the constraint 
operator is not symmetric then the group averaging method cannot be immediately 
applied. However, for the symmetric operator, to which we can apply the method, 
direct methods are not possible to apply and so group averaging becomes important. 
Expectedly, in the diffeomorphism invariant context it is true that the constraint 
algebra becomes Abelian which makes group averaging especially attractive (see [0]). 



4 Generator of the Wick rotation transform 

As explained in section 3, since classically the generator is given by C = (tt/2)K, 
almost no work is needed to do. We just notice that K = —{V/k, H e [1]} and now 
simply define 

K:=-^- 2 [V,H E [l]] a ndC:=^K . (4.1) 

Notice that V, H E have both dense domain and range in Cyl 3 („4/£). Therefore, K 
has also dense domain and range in Cyl 3 (*4/C?). Moreover, consider the case that 
we are using the symmetric version of H E (N). Then K is also symmetric and since 
fT^fl] is an imaginary- valued operator, K and V are both real- valued operators. 
We now choose complex conjugation as the conjugation in von Neumann's theorem 
(see [J!]]) and see that both V, K have then self-adjoint extensions. 
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Remark : 

Notice that the method of getting K from V, H E [1] simply through the commutator 
of these two operators is not possible to apply using any of the other operators 
corresponding to ii^fl] as defined so far in the literature : the classical identity 
K = — {V, H E [1]} holds only for H E given by ( j2.7|) , it does not hold either for 

\J det {q)H E |28|, |29| or for J ^Jdet(q)H E |14], |15| and thus for those approaches the 
operator corresponding to the generator of the Wick rotation transform is far from 
easy to define. 

For later purposes we wish to derive a more explicit expression for K. Let / be a 
function cylindrical with respect to a graph 7. Then we have 

Kf = -JS E \V^Hf\f = ~ £ [V v ,H E ]f=: £ K v f (4.2) 

where in the first step we again used the fact that the volume operator does not see 
the vertices V(7 U A (7)) — V(j) and in the second step we exploited that V v only 
acts on the edges incident at v so that the commutator with H E , which contains 
only holonomies of edges incident at v', vanishes if v' 7^ v. So, each term K v in the 
sum over vertices contains only the part V v of the full volume operator. 

5 The Wheeler-DeWitt constraint operator 

As outlined in section 3, we will now first derive a regulated operator corresponding 
to the expression tr([K a , K b ] [E a , E b ])/ y / det(g) which is consistently defined, whose 
action is diffeomorphism covariant and, in case we are dealing with a symmetric 
H E (N), is symmetric as well. In the symmetric case the idea is then to treat it 
as a perturbation of some self-adjoint extension of H E in the expression of H and 
to try to invoke the Kato-Rellich theorem to conclude that it is self-adjoint on the 
domain of H E (see |J). Finally we will prove that both versions of the operator are 
anomaly-free. 

5.1 Regularization 

We will not repeat all the arguments here as the derivation is completely analogous 
to the one for H E . 

We begin with the classical expression for the "kinetic term" 

r N 

T[N] := 8 / d 3 x-e abc tr({A a , K}{A b , K}{A C , V}) (5.1) 

and introduce the same triangulation T as used for H E to regulate the integral with 
the result 

T[N]t := ~ E WA))^tr(/^ ( ^^ 
AeT 

(5-2) 

Here, N is again the lapse function divided by k. Each term in the sum ( |5.2|) labelled 
by a tetrahedron A tends to 8/k 3 f A JVtr({A, K} A {A, K} A {A, V}) as we shrink 
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the tetrahedra to their basepoints, that is, we get the correct continuum limit. 
The quantization of this expression now merely consists in adapting T to a graph as 
outlined for the regularization of H E and in replacing Poisson brackets by commuta- 
tors and functions by operators. The result, when applied to a function cylindrical 
with respect to a graph 7, reads 

xtr(/i Si(A) [/i-| A) , K v )h s . {A) [h-\ A) , K v )h Ski A)[K^ A) ,V v ))f 

uev( 7 ) u(A)=« »ev(7) J ^y v > 

In order to arrive at this expression we use the same arguments as before to conclude 
that only tetrahedra intersecting a vertex of 7 give a contribution and we realize that 
all the commutators with holonomies along edges incident at v are non-vanishing 
only for the parts K v and V v of V and K respectively. To prove the latter we merely 
need to observe the following : 

Let s be any segment of an edge of 7 incident at a vertex v then g := [h s ,K]f = 
£t/ev(7) h a K v f - T,vev(~fU{s}) K v h s f = [h s , K v ]f because of the already familiar ar- 
gument that K v i only involves the volume operator at v' which commutes with 
h s unless v — v' and we also used s U 7 = 7. Next we have already seen that 
h := h aj (A)[hJ (&),V]g = fosj(A)[^V(A)' Vv]g which is a cylindrical function with re- 
spect to 7 U A (7). Now, since for v G 7 and v' G V(j U A (7)) — ^(7) we have 
trivially [h^ A) ,K v >}h = it follows that h 3i ( A )[hJ^, K\h = h s ^ A) [h^ A) , K v ]h as 
claimed. This furnishes the regularization part. 



5.2 Cylindrical Consistency 

Notice that in (|5.3| ) there appear factors of h a jM,A = A(7) which are holonomies 
along segments of edges of the original graph 7. Now, according to the definition 
of K one is supposed to adapt the triangulation associated with K (which is of 
course the same as the one associated with H E ) to the graph 7 U Sj(A(7)) when 
acting with K on h^^y^fy f° r instance. However, K is consistently defined up to 
a diffeomorphism and 7 U ei(A(j)) = 7 so that ( |5.3|) is cylindrically consistently 
defined up to a diffeomorphism, provided we introduce again the projectors pa- So 
we redefine 

t - % := Ta^~ (5.4) 

«(A)=A E ( v ) 

for the non-symmetric version of H(N) while for the symmetric version we refer to 
(3.5), (3.6) and (3.7) of 0. 



5.3 Diffeomorphism-covariance 

Diffeomorphism covariance is immediate because T is covariantly defined if and only 
H E is. There are no loop assignments different from the ones made for H E involved 
in the definition of T. 
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5.4 Anomaly-freeness 

Theorem 5.1 The complete (symmetric or non-symmetric) Lorentzian Wheeler- 
Dewitt operator 

H 7 [N] := f,[N] - H E [N] = £ N v [f v -H E ] (5.5) 

is anomaly-free. 
Proof : 

To see this we use the fact that neither H E nor T act at the additional vertices 
introduced by acting with H E , T (in the symmetric case this requires Lemmata 3.1 
and 3.2 and find that for a spin-network function / cylindrical with respect to 7 we 
have 

[H[M], H[N)]f — £ Mv "r~J?* Nv [kit, T v ,\ + [g* g*]} - fc, gfl]/ 

(5.6) 

that is, again the vertices of every V(7 U A (7)) — V(j) were irrelevant. In (|5.6|) it is 
understood that the operators T v , H E still depend on the graph on which the function 
depends that they are acting on in complete analogy with the considerations made 
in Q3.14[ ). Now we see that for v = v ' the factor consisting of the antisymmetric 
product of the lapse functions vanishes (which was crucial in replacing E(v') by 
E{y') = E(v',j)) while for v 7^ v ' the commutators vanishes when evaluated on a 
diffeomorphism invariant state as in after getting rid of some diffeomorphism 

operators, again because all operators involved in the commutators only contain the 
parts V V ,V V > of the volume operator which do not act on edges incident at different 
vertices. 

Let us be more specific. Notice that the action of the operator T^, v a vertex of 7 
such that v(A) = v , A = A(7), on / looks in full detail like this (c is a universal 
constant) 

^\ u i I J «(A( 7 ))= M (A'(7UA(7)))=u 

Xtr ^(4)[^(4). [^^^(^JAO^cao^Vo,^])]]^ x 

X K (A) K(A) ' eP9rtT (V, ( A ')^,(A) [^.KIPmA)^,'^)/' 

While this is a complicated expression, what counts is that the tetrahedron pro- 
jectors Pa'( 7 uA(7)) can be replaced by Pa'M, that is, a projector with respect to 
the original graph 7 and therefore also E(v,'j U A(7)) can be replaced by E(v,j). 
The same is true for the projectors that are involved in the operator T^, or H^, 

when applying it to T^f or H?f where v(A') = v' ^ v. The associated tetra- 
hedra involved depend on a nested sequence of graphs, that is, we have terms 
involving tetrahedra depending on 7 C 7 U A (7) C 7 U A (7) U A'() U A (7)) C 
7 U A (7) U A' (7 U A (7)) U A" (7 U A (7) U A' (7 U A (7))). However, just as in the case 
of computing the commutator between two Euclidean Hamiltonian constraints, we 
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see that we get a sum of terms each of which looks like U(tp)AB — U(ip')BA. This 
happens for each pair of values that the pair of involved operators (p^E(v),p^,E(v') 
may take because they act at different vertices and so it is irrelevant which acts first, 
they commute. This, together with the fact that the volume operator at one vertex 
does not act on the tetrahedra incident at another vertex shows that A, B commute 
upon removing the diffeomorphisms ip, ip'. 
This furnishes the proof. 
□ 

In the companion paper jy] we will compute the complete solution to both the 
non-symmetric Euclidean and Lorentzian as well as the Diffeomorphism constraint. 
Furthermore we define the symmetric version of the Lorentzian operator show that 
it is diffeomorphism-covariant and non-anomalous as well, however, we do not have 
the complete kernel in this case. Also we will report on the status of the Wick 
rotation transform in the light of these two articles. 

Finally, we will conclude with an interpretation of the results obtained, in particular, 
why the series consisting of these two papers was named "Quantum Spin Dynamics 
(QSD)". 
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